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Abstract. We consider the distribution of cycles in two models of random per- 
mutations, that are related to one another. In the first model, cycles receive a 
weight that depends on their length. The second model deals with permutations 
of points in the space and there is an additional weight that involves the length of 
permutation jumps. We prove the occurrence of infinite macroscopic cycles above 
a certain critical density. 
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1. Introduction 

Random permutations and their cycle structure have been studied for many decades, 
with a strongly increased activity in recent years. Apart from the rich mathematical struc- 
ture, this interest is justified by a wide range of applications, from Gromov-Witten theory 
[lOj to polynuclear growth |i7j and to mathematical biology [5j. Detailed properties have 
been established in the case of uniform permutations for the moments of the distribution 
of the n-th shortest (or longest) cycle [13], for the longest increasing subsequence [1], or for 
the convergence to equilibrium ^T^■ Physics and biology have suggested certain models 
with nonuniform permutations for which interesting results have been obtained |1H [6]. 

In the present work, we introduce a model for nonuniform permutations which is moti- 
vated by its connection to the theory of Bose- Einstein condensation |8l [151 13 12] • Mathe- 
matically, the distinguishing feature of our model, when compared to the works mentioned 
above, is that the measure on permutations possesses a spatial structure. More precisely, 
we consider pairs (a;, vr) with x G (A is a cubic box in M'^) and vr € Sn (the group of 
permutations of N elements). The weight of {x, vr) is given by the "Gibbs factor" e"^^"'''^-' 
with Hamiltonian of the form 

TV 

H{x, vr) = ^ ^{xi - x^(i)) + ^ aeriin). (1.1) 

i=l l^l 

We always assume that ^ is a function ^ M U {oo}, with J d'^x = 1. The cycle 

parameters ai, 02, • • • are some fixed numbers, typically but not necessarily positive, and 
r^(7r) is the number of cycles of length ^ in the permutation vr. The length of the cycle that 
contains the index i is the smallest integer £^ 1 such that 7r^(i) = i — this definition of the 
length involves the permutation but not the underlying spatial structure. Intuitively, the 
Gibbs factor restricts the permutations so each jump is local, i.e. the distances \xi — 
remain finite even for large systems. The main question deals with the lengths of the 
cycles, in the limit of infinite volumes. 
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When ^(x) = '^\x\^ + c and = for all we obtain the model of spatial random 
permutations that corresponds to the ideal Bose gas; in this case pc is the well-known crit- 
ical density for Bose-Einstein condensation for non-interacting particles. The occurrence 
of macroscopic cycles in the ideal Bose gas has been understood in [15l |T6]. The present 
setting with general functions ^ was considered in [2] . The latter article also introduces the 
cycle weight 02 as an approximation for the interactions between quantum particles, and 
the occurrence of macroscopic cycles was proved for large densities. The present article 
extends the results of [2] to more general cycle weights and to all densities larger than the 
critical density. 

The main result of this article deals with the occurrence of infinite permutation cycles. 
Let iVfc(7r) be the random variable that counts the number of points in cycles of length k 
(we have Ylk ^k{T^) = ^ for all vr), and let E\^N{Nk) be its expectation. We consider the 
thermodynamic limit lAjjA/" — > 00 with fixed density p = N/\A\ (|A| denotes the volume 
of A). Fatou's lemma implies that 

E|ir„^A..iAi(||)<p- (1-2) 



We prove in this article that the left side is strictly less than p if and only if the density 
is larger than a critic al d ensity pc G (0,cxd]. The precise formulation of this result can 



be found in Theorem 3.5 We need to restrict to certain functions ^ (namely, e ^ has 



positive Fourier transform) and small cycle weights, in the sense that — > as £ — > 00, 
faster than l/log£. Our results include an explicit formula for pc, cf. Eq. (3.8), and 
some characterization of the nature of infinite cycles — they are macroscopic and the 
distribution of a given macroscopic cycle is uniform in [0,p — pd- Points that are not in 
macroscopic cycles are shown to be necessarily in finite cycles, and their density is given by 
max(/9, Pc). To our knowledge, the presence of a density of points in finite cycles is specific 
to spatial models, and does not occur in the other known models of random permutations. 

The model considered here can be thought of as the "annealed" version of another 
"quenched" model where the space positions would be fixed. The latter model seems 
more attractive from a probability perspective, and it is actually discussed in [2 . There 
are two main advantages of the present annealed model: the availability of rigorous results 
about infinite cycles, and its closeness to the Feynman-Kac representation of the quantum 
Bose gas. 

The structure of this article is as follows. In Section [2] we introduce an auxiliary model 
of non-spatial permutations with cycle weights, which turns out to be closely related to 
the spatial one. It corresponds to taking x = (0,0,...). We only consider the case of 
small weights ag so that the typical nonspatial permutations are like those with uniform 
distribution. We will study more general weights and other behaviors in a subsequent 
article [3^. 

In Section [3] we introduce the model of spatial permutations with cycle weights. We 
discuss the existence of the infinite volume limits for thermodynamic potentials and the 
equivalence of ensembles. These notions belong to statistical mechanics rather than prob- 
ability theory, but we need these results when we consider the more relevant question, 
as far as probability theory is concerned, of the occurrence of infinite cycles. Our main 



result is Theorem |3.5| in Section 3.2 We find in particular that the cycle weights modify 
the critical density; so they do have an effect on spatial permutations, unlike what was 
observed in Section [H 
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Section |4] is devoted to the proofs of the results about the thermodynamic potentials. 
We adapt the classical methods of Fisher and Ruelle |12] to our context. In Section [s] 
we relate our spatial model to an equivalent model in the Fourier representation, which 
allows to prove our main theorem. The nonspatial model with cycle weights also plays a 
role here. We use techniques introduced for the ideal Bose gas in [4J and |16j . 



2. The simple model of random permutations with cycle weight 

We start with the study of random permutations of n elements with no spatial structure, 
but with cycle weights. The results of this section will be useful for the spatial model and 
this is our main motivation. But the present model has its own interest. We consider only 
small cycle weights here, but results for other regimes will be presented in a subsequent 
article |3|. 

2.1. Setting and properties. The sample space is 5„ and the probability of a permu- 
tation IT £ Sn is given by 

PnX^) = -r— rexpj- ^ a£r£(7r)| (2.1) 

with normalization 

hn = ^y^ e-^^°*''^W. (2.2) 

Here, ai,a2, ■ ■ ■ are fixed numbers and r^(7r) denotes the number of cycles of length i in 
the permutation tt. Notice the symmetry: since ^jjrj^n) = n for all vr, the probability 
Pn is invariant under the transformation 

aj I— > aj + cj, (2-3) 

for any constant c; the normalization satisfies /i„ i-^ e~'^"'/i„. 

Let Na^bi'^) = Yl'i=a ^^^(^) denote the number of indices that belong to cycles of length 
between a and b. Our main result deals with the asymptotic distribution of cycle lengths. 

Theorem 2.1. If > i 1|1 ~ e^"* | < oo, we have for any ^ s ^ 1 

lim l^EniNi^sn) = s. 
n~*oo 

In essence, the hypothesis of the theorem requires that ai ^ a bit faster than 1 / log i. 
Theorem |2.1| implies that almost all indices belong to cycles whose length is a positive 
fraction of n. It can be shown that the number of cycles is of order log n. The claim is 
easy to get in the case of uniform random permutations (a^ = 0), but the extension to 



even small weights requires some efforts. The key to the proof of Theorem 2.1 which is 
given in the next subsection, is the following relation. 



Lemma 2.2. 

b 



For the proof of this lemma, we just remark that En{Nafi) = npn{^i G [ffl)&])) with 
ii = ^i(7r) the length of the cycle that contains 1. Summing over all possible values j of 
ii, and observing that there are possible cycles, we get the relation above. 
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2.2. Properties of the normalization In view of Lemma 2.2 it is clear tliat we 
need to gatlier some information on hn- We start witli a few exact relations (Proposition 



2.3) and we then obtain estimates (Proposition 2.4). 



Proposition 2.3. The hn 's satisfy the following properties 
(a) A recursion formula: 



hn = -y^e hn-e, ho = 1. 

(b) An explicit formula: 



n 



n ^ k g— Cff. 



A;! 

k=l ii,...,£k ^ 1 «=1 
£i+...+4=n 



(c) The increments satisfy 



/i„-/i„_i=^- Yi n — i 

k=l ' h,-/k ^ 1 i=l * 



(d) Another formula for hn: 



1 ^ A. e-"^« - 1 



"-E^ En,- 

fc=0 ^ 1 i=l 

li + ...+lk ^ n 

(e) If ^ ^ and 7 > 0, we /ia?;e 

E e-^"/in = exp J] ^ . 

The formula (b) shows that /i„ is decreasing with respect to (a^). The formula (d) 
shows that /i„ is increasing with respect to n if ^ 0. Now we cannot resist but ask the 
reader to consider the following expression: 



"1 1 
k=i ' ^1 • • ■ 

How does it behave for large n? The answer is surprisingly simple and is given by Propo- 



sition 



2.3 (b). 



Proof. The recursion formula is obtained from Lemma 2.2 by noting that En{Ni^n) = 
For the claim (b), it is useful to define bj = e~°'^+^ , j ^ 0. The recursion formula can be 
written as 

n 

(n + l)/ln+l = ^i^n-j- (2.4) 

Thus the series ((n + \)hn+i) is equal to the convolution of the series (6„) and {hn). We 
introduce the generating functions 

Gh{s) = Y hns"", Gb{s) = Y bus"". (2.5) 

n>0 n > 
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The generating function for the series ((n + is G^(s). By the properties of convo- 

lutions, we have 

G'^is) = Gb{s)GH{s), Gh{0) = 1. (2.6) 

The solution is 

Gh{s) = exp r Gb{t)dt = exp V (2.7) 

Expanding the exponential and rearranging the terms allows to find an expression for each 
coefficient hn- This gives Proposition 2.3 (b). 

The generating function for the increments, Gsh, satisfies 

Gsh{s)= ^{K-K-i)s'' = {l-s)Gh{s) (2.8) 

n > 



(with h-i = 0). Using (2.6), we get a differential equation for Gsh, namely 

G'sds) = Gsh{s)[G,{s) - j^] (2.9) 

The expression in the bracket is equal to the generating function Gb-i of the series (6„ — 1) • 
Solving the differential equation, we get 

Gshis) = exp Gb-i{t)dt = exp V (2.10) 
Jo n 



n > 1 



Again expanding the exponential and matching the coefficients, we get Proposition 2.3 



(c). The formula (d) follows from (c) since /i„ = Yll=o{^i ~ hi-i)- Finally, (e) follows 
directly from (2.7). □ 

We now collect a few estimates for 

Proposition 2.4. 

(a) // J2e ^ 1 z|l ~ I < oo, we have ^„ \ hn — hn-i\ < oo, and 

-„ • 

1 „— -nai 



(b) K ^ liK-ie "1 ; iterating, K ^ 

(c) ^e-"". 

(d) If {at) is subadditive (i.e. ctY^ ^. ^ 'Yli^ti)' ^ e" 

(e) If (ae) is superadditive, hn ^ e~"" . 



Proof. The claim (a) is an immediate consequence of Proposition 2.3 (d) and the domi- 
nated convergence theorem: In the limit n — > oo, the constraint ii + . . . ik ^ n vanishes 
and the corresponding expression factorizes. (b) and (c) follow from the recursion formula 
of Proposition |2.3| (a), keeping only the term £ = 1 for (b), and the term £ = n for (c). 



For (d) and (e), we write the formula of Proposition 2.3 (b) as 



^ 1 ^ e-E?-"', 



kl ^ £\...£k 

k=\ ^ 
^l+■■■+4='^ 



We replace ^ at,, by a„, getting an upper bound for hn if (a^) is subadditive, and a lower 
bound if [ai) is superadditive. □ 
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Proof of Theorem 2.1 We split the expression of Lemma 2.2 into 



i=i i=i 

We suppose that sn is an integer; it is easy to adapt the proof otherwise. Since — > hoc, 
for any £ > there exists such that if n — j > n^, 

1-e <l + e. (2.13) 



The last term of (2.12) is then less than e if s < 1 and n large enough (the case s = 1 is 
trivial). For any 6 > 0, there exists ns such that X]j>ni jl^ ~ e""-' | < 6. The sum over 



the first ng terms in the middle term of (2.12 1 is then less than (1 + s)^ j|l — e 



and it vanishes in the limit n — > oo; the sum over the remaining n — ns terms is less than 
(l + e)(5. □ 

We conclude this section with a corollary which follows immediately from Proposition 



2.4 and which will be very useful later. 



Corollary 2.5. Let 

h 

T-> ™' 

B = sup — — . 

m,n ilfi 

Then if > ||1 — e~"^ | < oo, we have < B < oo. 

3. The spatial model of random permutations 

We now introduce the spatial structure and consider the model described in the intro- 
duction. Let A be an open bounded subset of M'^, and let be the number of "particles" 
of the system. The state space of our model is 

Qj^^^ = A^xSn (3.1) 

with Sjy the group of permutations of N elements. Qa,n is equipped with the product 
of the Borel a-algebra on A^, and the discrete c-algebra on S]\f. Let {x,7r) G ^a,n with 



X = {xi, . . . ,xn)', our Hamiltonian is given by ( 1.1 ). 



3.1. The thermodynamic potentials. We consider the canonical ensemble where the 
particle density is fixed, and the grand-canonical ensemble where the chemical potential 
is fixed. The canonical partition function of this model is 

nA,iV) = ^ [ dxY^ e-^(-'-). (3.2) 



iV! 



The division by A^! guarantees that the partition function scales like the exponential of the 
volume. The definition makes sense for integer N; it is convenient to extend the partition 
function to noninteger A^, e.g. by linear interpolation. The grand-canonical partition 
function is 

Z(A,/x)= e''^y(A,Ar). (3.3) 

N >0 
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The parameter // is called the chemical potential. Then we can define the two relevant 
thermodynamic potentials, the free energy and the pressure: 

qA{p) = -T^AogY{A,\A\p), 

I I (3.4) 
Pa(m) = |^logZ(A,^). 

The parameter p is the density. 

We define the "dispersion relation" e{k) by 

^('^^ = / e-^'^''^^ e-«(^) dx. (3.5) 



e 

For now £{k) can be complex; e(0) = and Ree(A;) ^ a\k\'^ for small k. The most relevant 
case is the Gaussian, e~'>^^^ = (47r/3)~''/^ e"'^' . This corresponds to the ideal Bose gas 
and £{k) = 47r^/3|/cp. We will suppose in the next subsection that e(fe) is real. 
Note also that the transformation aj ^ aj + cj translates into 

9a(p) ^ QAip) + cp; 
Vh\N ^Pa{p-c). 

Next we recall the notion of Fisher convergence |I2^. A sequence (A„) of domains in M'^ 
converges to M'' in the sense of Fisher if 

• lim^^oo |An| = oo. 

• As e 0, 

f^ediamAn „ 
sup r-r— ^ > 0, 

where drA = {x G M'^ : dist {x, dA) ^ r}. 

This notion is very general. If A is bounded with piecewise smooth boundary, then the 
scaled domains A„ = {nx : x € A} form a Fisher sequence. 

For the following three theorems, we always suppose that a^/i converges as £ — > oo. 
Because of the symmetry (3.6) we can choose the limit, so we suppose that 

lim ^ = 0. (3.7) 

The first result is about the infinite volume pressure, that is given by an exact expression. 

Theorem 3.1. For any G M \ {0}, and any sequence A„ that converges to M"' in the 
sense of Fisher, we have 

Notice that p{p) is finite and analytic for p < 0, and that p{p) = oo for p > 0. 

Next, we have the existence of the thermodynamic limit for the free energy. 

Theorem 3.2. There exists a convex function q{p) such that, for any sequence (A„) of 
domains converging to (Fisher), and any sequence (pn) of numbers converging to p ^ 0, 
we have 

lim q\„{pn) = q{p)- 
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Then q is continuous. It is a standard exercise in analysis to show that the above prop- 
erty is equivalent to uniform convergence of qn to q on compact intervals. Finally, pressure 
and free energy are related by Legendre transforms, a property known in statistical me- 
chanics as "equivalence of ensembles" . 

Theorem 3.3. The infinite volume pressure and free energy are related as follows: 

q{p) = sup[pfi-p{n)], 

p{p) = sup[p/i-g(p)]. 
p 

And q{p) is analytic except at the critical density 



E 

n > 1 



-ne{k) 



(3.8) 



One can check that the critical density is real. It is always finite in dimensions d ^ 3; 
but it may be infinite in d = 1, 2, in which case q is real analytic for all p in [0, oo). 








(a) 



(b) 



(c) 



Figure 1. Qualitative graphs of the pressure p{p) and its Legendre trans- 
form the free energy q{p). The critical density pc is infinite in (a) and (b), 
and it is finite in (c). 



So far we have considered free boundary conditions. For the proofs of our results on 
the probability model below, we will need versions of the above theorems with periodized 
boundary conditions. To be precise, let us consider a sequence (A„) of d-dimensional boxes 
with side-length L„, with L„ — > oo. We define i^A as in (3.9) below, and let q^^{p) and 
p^^'^(^) be the corresponding free energy and pressure. As is usual in statistical mechanics, 
a change in boundary conditions brings a correction to thermodynamic potentials of the 
kind g'^*"^ — ^a ~ with \dA\ a measure of the boundary of A; same for the pressure. 
The next theorem is less sharp but it is enough for our purpose. 
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Theorem 3.4. For any p ^ and any G M \ {0}, 

n—^oo " 

lim PaJ(^) =p(^). 
In both cases, convergence is uniform on compact intervals. 



While Theorems 3.1f 3.3 are proved in Section [4j Theorem 3.4 needs the notation of 



Section [5] and it is proved in the appendix. 

3.2. The probability model. We will now study our model from a probabilistic point 
of view, proving the occurrence of infinite cycles above the critical density. As is often the 
case, we will need more stringent conditions than we did for studying the thermodynamic 
potentials. 

We also need to slightly modify the Hamiltonian. Let A be a d-dimensional cubic box 
with side length L and volume V = L'^. Define through 



The important point is that e^^* has positive, A-independent Fourier transform e"^'-'^-* al- 
ready in finite volume; t his h elps us to relate our model with a probability model on Fourier 
modes. See Proposition 5.1 in Section [sj In particular, we note that e"^''^^) d'^x = 1. 



If e~^ has compact support and if L is larger than the diameter of the support, this 
"periodized" setting corresponds to usual periodic boundary conditions. 
We now define 



N 



HA{x,7r) = "^Uixi + ^ a^r^(vr). (3.9) 



i=l i>i 



and introduce a probability measure on such that a random variable 6 : rJA,7V ^ 

has expectation 

Note that y(A, N) is the partition function with periodized boundary conditions, a fact 
that we suppress from the notation. 

Next, let ^j(vr) = 1,2,... denote the length of the permutation cycle of vr that contains 
the index i. It is convenient to consider the density of points that belong to cycles of 
certain lengths. Precisely, let 

Qa,,{^) = ^#{^ = 1, 2, • • • : a ^ ^(vr) ^h} = (3.11) 

Theorem 3.5. Assume that e~^ is continuous, that it has positive Fourier transform 
(i.e. £{k) is real), and that > ^ ^ < oo. We also suppose that pc < oo. Let rj be any 
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function such that rj(y) — > oo and 'i](y)/V ^0 asV^ oo. Then for all s ^ 0, 



lim EA^pviOi^rjiV)) , .r ^ 

V^co \^pc tf Pc 



(microscopic cycles) 
(mesoscopic cycles) 



lim EA,pv{Qn{v),sv) 

V—*oo " ' 



ifp^pc, 

s if ^ s ^ p — pc (macroscopic cycles) 

P-Pc ifO^p-pc^s. 



Recall the expression (3.8) of the critical density. It is easily seen that 

dk 

g£(fc) _ I 



Pc 



oo <^ 



OO 



dk 
\<i e(^) 



oo. 



(3.12) 



The same theorem is stated in ||2J but with ae = 0; in this case it holds also when pc = oo. 
The proof of Theorem 3.5 , which can be found at the end of Section|5] requires Proposition 



5.4, which states that certain occupation numbers of Fourier modes are typical. But we 



can prove one of the claims of Proposition 5.4 only if < oo. There is little doubt that 
all properties also hold true when the critical density is infinite. 

4. Thermodynamic limits and equivalence of ensembles 

In this section we consider the thermodynamic limits of the pressure and of the free 
energy, and we prove Theorems 3J_-3^, We adapt the methods devised in the 1960's by 
Fisher and Ruelle for classical particle systems. See [12] for references. 



Proof of Theorem 3.1 The case 02 7^ 0, a„ = for any n ^ 2, was treated in [2]. The 
general case is similar. The key observation is that everything factorizes according to the 
permutation cycles. 

The grand-canonical partition function reads 

Z{A,f^)= nil dxi...dx„e"^"=i«(^'"-^'+i) e-"")'"^"^ (4.1) 

where r„(7r) is the number of cycles of length n in the permutation vr. We also supposed 
that Xn+i = xi. The number of permutations of elements with r„ cycles of length n, 
n ^ 1, is equal to 

iV! / JJ n''"r„!. 



n > 1 



Then 



^(^'^)= n{E 

n ^ 1 r ^ 



1 / & 
J' I 



pn-an 



n 



dxi . . . dXn Y\ 6 



i=l 



exp< 



^pn-an 



n 



n 

dxi...da;„ JJ e-«(^'-^'+i) }, 



(4.2) 



i=l 



and thus 



Pa(m) 



lAI 



E 

n > 1 



n 



dxi . . . dxn e~ 



-((Xi-Xi + l) 



(4.3) 



i=l 
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Now, with xi = in the second term, 

„ n 

-S,{Xi-Xi + l) 



.\r [ dxi...dx,ne-«(^»-^»+iu / dx2...dx„n 

1^1 J A" fJl JMd{"-i) fJl 

\0) (4-4) 



The inequahty becomes an identity in the hmit A ^ R'^, and therefore we get Theorem 



3.1 by dominated convergence if /x < 0. Since e(0) = and e{k) > a\k\'^, the integrals of 
g-?ie(fc) around A; = decay as an inverse power of n, and the integral for large k decays 
exponentially. Then if ^ > 0, we have PA^fJ-) = oo for any A large enough. □ 

Let F{A, N) denote the macroscopic free energy, 

F(A,iV) = -logy(A,A^). (4.5) 

Lemma 4.1. Let Ai,A2 be disjoint open bounded subsets ofW^. Then for any integers 
Ni,N2, 

F(Ai U A2, iVi + N2) ^ i^(Ai, iVi) + F{A2, N2). 

Among the many useful consequences of this subadditive property, we get the upper 
bound: 

TV 

F{A,N)^ ^F(Ai,l) =iV(e(0) + ai + log/5), (4.6) 

i=l 

where the Aj's partition A in subdomains, each of volume |A|/A^ = 1/p. 



Proof. We show that Y{Ai U A2, iVi + iV2) ^ ^(Ai, iVi)y(A2, iV2). From the definition, 

y(AiuA2,iVi + iV2)= i/vv / '^y H 

(iVi + N2)l y(A,uA,)-i y(A,uA,)-2 ^^^^^^^ 



(4.7) 

We get an upper bound by restricting the integrals so that exactly A^i points fall in Ai 
and points in A2. Rearranging the integrals, we obtain 

y(A,uA2,A^, + iV2)^^/^^^d./^^^dy e--(--). (4.8) 

We restrict the sum to permutations of the kind tt = tti x 112, with tti a permutation of 
the first A^i elements and tt2 a permutation of the last N2 elements. Then 

y(A, U A2, iV, + iV2) ^ ^ 1^^^ dy e-^(— -) . (4.9) 

1 2 7rie5jvj^ 

•'r2 6'Sjv2 

Finally, we observe that 

H{x X y, TTi X 112) = H{x, TTi) + H{y, 112). (4.10) 



Integrals in (4.9) factorize, yielding the product of partition functions. □ 
Next we identify the free energy by considering a special sequence of increasing domains. 
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Lemma 4.2. Let Cn be the cube of size 2" centered at the origin, and define 



q{p)= lim 2-'^"F(a,L2>J). 



Then 



(a) the above limit exists indeed. 

(b) q{p) is convex. 

Proof. The existence of the hmit fohows from a standard subadditive argument, where we 
show that 

Hmsup2-'^"F(Cn, [2"/)J) ^ hminf 2-'^"F(C„, [2>J). (4.11) 

n — ^oo ^ 

Fix k; we consider n > 2k. The cube C„ can be partitioned into 2'^^'^^^^ cubes C^. We 
put [2'^^p\ points in the first 2'^("-'=) - 2'^("-2fe) ^^^^gg^ ^^le remaining Nnk points in 
the remaining 2'^^"^^'^) cubes. By subadditivity, and the upper bound (4.6), we get 



F(C7„,L2^>J) ^ (2^(^ 



n— fc) nd{n—2k)\ 



)F{Ck, L2'^Vj) + A^nfc(e(0)+ai + log 



Nnk 



■ (4.12) 
(4.13) 



2fi(n-fc)_ 

The number A^^^ is not too big: 

Nnk = L2'^"/5j - (2'^(''-'=) - 2'^('^-2'=))[2'^VJ ^ 2'^("-'=)(l + p). 
Then 

hmsup2-'^-F(C„, L2>J) ^ 2-'^'=(l-2-*)F(C,, [2*/,J )+2-^^(l+p)(C(0)+ai+log(l+/5)). 

n— »oo 

(4.14) 

Taking the liminf A; ^ cxd in the right side, we obtain (4.11). 



We now prove that q{p) satisfies a certain form a continuity, see (4.20) below. Let e > 0; 
we study 

g(p + e) =lim2-'^'^F(C„, [2'^'' {p + e) \) . (4.15) 

n 

Fix k. There are 2'^("-'') cubes Ck in C„. Put [2'"'pJ points in the first [2'^("-^)(l - e)\ 
cubes, and the remaining Nnk points in the remaining domain Dnk- By subadditivity and 
the upper bound ( |4.6| ), 

F{Cn, L2""(p + e)J) ^ l2''(^~'Hl -e)\F{Ck, L2*pJ) + iV„fc(C(0) + «i + log ^). (4.16) 

We can estimate the last term: 

Nnk = L2'^"(p + s)\ - [2^("-^)(l - e)J [2"^^] ^ 2'^"+ie + 2'^("-^) + 2*p, 



\D 



nk I 



tdn. 



(4.17) 



Then 



A^nfe(^(0) + ai+log|-^) ^2'^"(2e + 2-^ + 2-'^("-^)/,)(C(0) + «i + log(2+2^ + 2zf^ 

(4.18) 



We substitute these bounds into (4.16), divide by 2 , and let n — > oo. We get 

q{p + e)^ 2-^'=(l - e)F{Ck, L2*pJ ) + {2e + 2'''){m + oi + log(2 + 2^)). (4.19) 



This is true for any fc. As A; ^ oo, 

q{p + e) ^ (l-e)g(p) + 2e(e(0) + ai+log2). 



(4.20) 
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Finally, we prove a weak version of convexity. Combined with (4.20), it implies that 
q{p) is convex indeed. For any pi and p2 such that 2'^"~^/9j is integer for n large enough. 



'\pi + \P2) = Iim2-'^"F(C7„,2'^"-Vi + 2'^"-V2 

n 

^ lim2-°'"r2^-iF(C„_i,2'^("-iVi) + 2'^"^^(Cn-i,2'^(""^) 

= \<l{pl) + \(l{P2)- 



P2) 



(4.21) 



Indeed, we put a density pi of points in half the cubes Cn-i, and a density p2 in the other 
half, and we used subadditivity. □ 



We can now prove Theorem 3.2 about the convergence of the thermodynamic limit for 
the free energy. 



Proof of Theorem 3.2. First, we show that 

1 



limsup ---F(A„, |A„|p„) ^ q{p). 



(4.22) 



Choose a cube C, and a number p' > p. Given n, let M„ be the largest integer such that 
Mn\C\p' ^ |A„|/3„. We pave M'^ with translates of the cube C. The volume of cubes inside 
A„ is at least |A„| — |5diamC-'^n|- Thus for n large enough, the number of cubes inside A^j 
is at least 



I A. 



l^diamC-^nl ^ l-^n | Pn ^ 



\C\ - \C\ p' 

By subadditivity of the free energy, 

F(A„, \K\Pn) ^ MnF{C, [\C\p\ ) + F{D, N) 



with 
Now 

Then 



F{D,N) ^ iV(ai +log ^ 



\D\'- 



N = \An\pn - Mn[\C\p\ ^ |A„|p„(l - ^ + ^) + \C\p 

|Z)| = |A„|-M„|C|^|A„|(1-^). 



(4.23) 

(4.24) 
(4.25) 

(4.26) 



Pn{l 



1 N , lev 

\C\p') + |A„| 



ai + log 



Pn(l- -7+177177) + 



\c\ 



p'-^|C|p'^^|A„| 



(4.27) 



Letting n — > oo, the expression simplifies a bit: 
limsupp^F(A„,|A„|/,„) ^ ^^F(C, L|C|pJ) + p(l 



p' \cv 



i-7)(ai + log 



\C\P' 



p(i- 



'p'~^\C\p'' 

1-4 
p 



(4.28) 

We can consider the cubes of Lemma 4.2 and take the limit k oo. Then we let 
p' p, and we get (4.22). 

We complete the proof by showing the complementary lower bound, namely 

liminf n^F(A„, |A„|/j„) ^ q{p). (4.29) 

n— >oo l-"^"! 
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A consequence of the limit in the sense of Fisher is that, given A„ ^ M'^, there exists 
r/ > such that each A„ contains a translate of a cube of size 77diamA„. Given n, let k be 
such that 2^^ ^ diamA„ > 2^"-*^. Then A„ is contained in a translate of C^, and 

1^1^^ 2-'^'=(?7diamA„)'^ ^ 2''^'=(r?2'=^i)'^ = {\t]Y . (4.30) 
By subadditivity, we have 

F(Cfc, [2^^^]) ^ i^(An, |A„|p„) \ A„, [2^^^] - |A„|p„). (4.31) 

As n ^ oo, we also have k oo. Then, taking the liminf of the expression above, we get 

liminf ^F(A„,|A„|p„) ^ liminf [|^2-'^^F(Cfe, yi"^ p\) 



-|A„| 



F(Cfe\A„, L2''VJ -|A„|/9„) . (4.32) 



|An| 



-|An| 



Notice that p\ - |A„|p„ = |Cfc \ A„|p^ with 



Pn 



L2''fcpJ-|An|p„ 
|Cfc\An| 



P 



(4.33) 



as n ^ oo. This uses (4.30). Using the result for the limsup, Eq. (4.22), we find that the 
right side of (4.32) is larger than (?(p). □ 



There remains to prove Theorem |3.3[ It is actually enough to show that p is the 
Legendre transform of g — since q is convex, it is necessarily the Legendre transform of p. 



Proof of Theorem 3.3. From (3.3) and ( |3.4[ ), we have 
Pk{p) 



(4.34) 



^log ^ exp||A|[/)/i-gA(p)]|. 

' ' peN/|A| 

By restricting the sum over a single /o, we get 

Pa{p) ^ Pl^- qxip)- (4.35) 

This holds for any A and any p such that |A|p is integer. It follows that for any p > 0, we 
can use a suitable Fisher sequence (such as cubic boxes of size L = p-^/'^k, ken) so as 
to obtain 

p{p) ^ pp-q{p). 

This inequality also holds when taking the supremum over p in the right side. 



(4.36) 



The upper bound requires a bit more work. By (4.35), we have 

- QAip) = \PP+ \PP - QAip) ^ 5PAt + PA(f )• 
Let ^ be a number that is independent of A, to be determined later. Then 



peN/|A| 



(4.37) 
(4.38) 



peN/|A| 
p>A 



The latter sum is equal to el^lPA(f ) (i _ )-i. 

It follows from Theorem 3.2 that qA{p) converges to q{p) uniformly on compact intervals. 
Thus for any e > 0, there exists A large enough so that 



pp - qAip) < p*p - q{p*) + e, 



(4.39) 



SPATIAL RANDOM PERMUTATIONS WITH SMALL CYCLE WEIGHTS 15 
for any < p < A. Then 

PAifi) ^ p*fi-q{p*)+e + j^^ log{^|A| + (1 - e^MlPAi^)-p'f^+<i{p*)-^HA^] |_ ^^ ^q^ 

We can choose A large enough such that the exponent is negative. Letting A M'^, we 
see that p{p) ^ snpp[pp — q{p)] indeed. □ 

5. The Fourier model of random permutations 



We consider now the setting of Section [3.2[ Thus A is a cube of side length L and volume 



V = L"^. Let A* = f Z*^ denote the dual space. The goal of this section is to describe a 



probability model on = (A*)^ x S^, whose marginal distribution on permutations 

coincide with the model of spatial permutations. The method is inspired by jl6], who 
studied the nature of cycles in the ideal Bose gas. The beginning of this section is almost 
identical to [2j, the difference being that ai ^ here. 

For a bounded domain A the new probability space is discrete. We give ourselves a 
positive continuous function e{k), k G M°'. We suppose that e(0) = 0, and that e{k) ^ a|A;p 
for small k. The probability of (fe, vr) G Q\ ^ is defined by 

|,A,^(fc,^)= W)^! (5.1) 
[0 otherwise. 

Here, the Hamiltonian H is given by 

N 

a£ri{TT). (5.2) 

1=1 e^i 

This model offers an alternate representation to the model of spatial permutations, as 
far as the permutations are concerned, thanks to the following relation. Let 6ij denote 
Kronecker's delta symbol. 



Proposition 5.1. Suppose that e ^ = e ^, and define Hj^ through (3.9). Then, for any 
permutation vr, 

■^^"^ fce(A')^ i=i 

In particular, it follows that y(A, N) = y(A, iV), and E\^m{J^^ x {vr}) = ps^^n{^Xn-> ^) 
for all permutations vr. 

Proof. We use the following identity, which follows almost directly from well-known rela- 
tions between Fourier transform and convolution: 

. N N 

See [2], Corollary 5.3, for more details. Notice that the case = 1 reduces to Poisson 
summation formula, and it holds true because e~^ is continuous. Multiplying both sides 
by e~S"*^**^'^) , one gets the result. □ 
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Next, we introduce occupation numbers. Let A/a be the set of sequences n = (11^) of 
integers indexed by A; G A*, and let Ma,n the set of occupation numbers with total number 
N: 

Ma,n = {n £ Ma : = nY (5.4) 

keA* 

To each k E (A*)^ corresponds an element n G Ma,n, with counting the number of 
indices i such that fcj = k. Thus we can view n as a subset of (A*)^. The probability K^ 
yields a probability on occupation numbers. Indeed, summing over permutations and over 
compatible vectors k, we have 

PA,N{n) = T77^ n ^""'^"^'^ (5.5) 



with hn defined in (2.2). 

We will obtain some properties on the probability pA,Ar(^) below. But we first relate 
this probability with the lengths of permutation cycles. 

Proposition 5.2. 

Proof. Let pA,N{k) = YlTTPA,N{k,7r). Since g^f, depends only on permutations, we have 
EA,N{Qab) = Y PA,N{k) Q^k{ir)pA,N{Mk). (5.6) 

fcG(A*)^ 7re5jv 

Actually, Qabi'^) depends only on the conjugacy class of vr. In other words, we have, for 
any vr, a S Sn, 

Qabi^^'^T^cr) = Qabi-^)- (5-7) 

It follows that 

Y Qab{T^)PA,N{Mk) = Y ^ab{'^)PA,N{'^\cr{k)). (5.8) 

Summing first over occupation numbers and then over compatible fc's, we get 

EA,N{Qab) = Y PA,N{n) Y QabMPA,N{Tr\k). (5.9) 

Here, k is any vector that is compatible with n. 

A permutation vr such that 7r(fc) = k (that is, = ki for all i) can be decomposed 
into permutations {TTk)k£A*, where ttj. is a permutation of the indices i such that ki = k. 
Notice that 

A^afe(vr) = Yl NabiTTk), r^vr) = Y ^^(^fc)- (5-10) 
keA* fceA* 

Then 

PAMk,n) = y^j^^_U e-S^-^'-^W (5.11) 

if 7r(fe) = k; it is otherwise. Also, 
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Then 

PA,iv(7r = — 77— = M — (5.13 

if 7r(fc) = fc, and otherwise. Using ( |5.10[ ), 

^ ^ Q—^lOnri(TTy) 

EA,N{Qab) = Y Yl P^'N{n) Nab{7^k) Yl 

neA^,iv (TTfcS^nJfceA* fc'eA* 



(5.14) 



□ 

In the hght of Proposition |5.2[ we can now focus on the quantity PA,Af(^)- Namely, 
our results from Section |2] imply that macroscopic cycles appear if and only if at least one 
mode is macroscopically occupied, i.e. iff p\^N{nk ^ sN) > uniformly in G N and A 
such that N = pA. 

We prove now that macroscopic occupation can occur only for k = and that it occurs 
if and only if p is above the critical density pc defined in (3.8 ). The first step is a result that 
gives detailed information about the limiting distribution of the random variable no/V. 



Theorem 5.3. Let po = ™ax(0,p — pc), with pc the critical density defined in (3.8). Then 
for all A ^ 0, we have 

lim ^A.pv(e^"«/^) = e^'^o. 

V^oo 

Proof. Our proof is based on the work of Buffet and Pule [3^ for the ideal Bose gas, see 
also We need to modify it due to the presence of cycle weights. We define 

Then Y{A, N,0) = Y{A,N). Recall that hoo = hm„^oo hn is given in Proposition 2.3 (d) 
and that < hoc < 00. Then by (5.5), 

P^Mno ^ J) = yr^ E 11 ^"""^'^ 

1 ^ /W. Yl e-"^^W/i„, (5-16) 



Y(A, N) 4^ 

Y{A,N-j,j) 
Y{A,N) 



for all J ^ 0. Using 



PA,N{no = j) = PA,N{no ^ j) - pA,Af(^^o ^ i + 1), (5.17) 
the change of summation index j ^ N — j gives 
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Here we used the convention ^(A, — 1, + 1) = and the fact that PK,N{nk ^ + 1) = 0. 
We now fix p ^ and put u = \/V . As a first step we show that 

hm SA,py(e^"°/^ ) = 1 for all p ^ pc- (5.19) 

Above, we wrote pV instead of [p^J , and we will continue to do so in order to simplify 
the notation. To prove (5.19), we get from (5.18) that 

Ek,pv{q V ) - q V = 



Y{K,pV) 
Y{A,pV) 



J^e-^ (Y{A,j,pV-j)-Y{A,j-l,pV-j + l)) - e" 

j=0 
pV 



py pV V 

e-^ y(A,i, pV - j) - E ^(^'■?' - ) 

0=0 j=0 ^ 



(5.20) 



(1 



_ A 

e V 



pV 



Y,e-vY{K,j,pV-j) 

j=0 



We need to show that (5.20) converges to zero as 1^ — > cxd. Recall the constant B defined 

we have 

e vY{A,j) 



in Corollary [2^5} Since y(A,A^,j) ^ BY{A,N), we have 

, pAp 



XriQ 

EA,pv{e^ 



e~v < B 



Y{K.pV)V 



j=0 



Be 



Xp 



V 



{l-e)pV 
j=0 



Aj/V Q-ViqA(j/V)-qAip)) 



epV 

+ E 

3=0 



'\{pV-j)/V 



Y{A,pV-j] 
Y{A,pV) 



(5.21) 



for any e > 0. Above, recall that (7a(j/^) is the finite volume free energy given by (3.4). It 

that q\ converges uniformly on compact intervals to the convex 
q{p) is strictly decreasing for p < pc- For each e > 



follows from Theorem 3.4 



function q. By Theorem 



3.3 



there is 65 > such that qA{j /V) — Qa{p) > for all V large enough, and all j ^ {l — e)pV. 
So the first term in the bracket above is bounded by (1 — e)pV e~''^^ and thus converges 
to zero as y ^ 00. For the second term, we claim that Y(A,pV — j)/Y(A, pV) ^ B for 
all j. This is proved by putting the extra j particles into the zero mode A; = 0, or more 
formally through 



y(A, AT) ^ E n ^ ^ - 



(5.22) 



no ^ j 



Thus the second term, along with the prefactor e^P \/V, is bounded by B^ Xe^P pe. As e 
is arbitrarily small, we obtain (5.19). 



We now turn to the case p > pc- We define the atomic measure 

00 

^^^,P = Ca E(y(A,i, pV - j) - Y{A,j -i,pV-j + l))5,/v 

j=0 



on 



with 



(5.23) 



(5.24) 
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here, we set 

AC\ = {n G AAa : no = 0}. (5.25) 

Later, we will also use the notation Ma.n = {n £ A/a, at : tiq = 0}. As we will see below, 
C\ is a correct normalisation so that /iA.p converges in the limit y ^ oo; it does not 
depend on p. 

We rewrite (5.18) using fiA,p, which gives 

J l[o,p](2;)MA,p(dx) 

The strategy is to study the Laplace transform of pA.p- It will be possible to take the limit 
V — > oo. Putting the limiting measure in the right side of (5.26), we will get the infinite 
volume limit of the left side. An advantage of this strategy is that several convergence 
issues are handled using standard theorems of analysis. The Laplace transform of /iA,p is 
given by 

oo 



poo _ _ 

/ e-"^ /XA,p(dx) = Ca(1 - e- ) V e-"^/^ Y{A, j, pV - j) 

•^0 j=0 

^ . (5.27) 

= C'A(l-e-VV) ^ (n e-(^«+VVK/,„ ^M^l^ 



neAfA ^feeA* ^ ■ 

with 

[hoo if EfceA* '^fc > ^• 

At this point, the idea in [4J and [2j was to factor out the contribution of the zero Fourier 
mode. This is not possible here because the factor /i(n, pV) couples the modes. This 
difficulty can be circumvented by introducing 

oo 

i^A = Ca J^(r(A,i, -1) - Y{A,j - 1, -1))<5,/^ (5.29) 

j=0 

and by writing 



oo />oo />oo 



e 



PA^dx) = / e-^^ i.A(dx) + / e-^^ (/XA,p(dx) - Mdx)) . (5.30) 



We now prove that the second term vanishes as F — > oo. We have 
/■oo ^ °o ^ 

/ e-^^(/iA,p(dx)-i.A(dx)) =Ca(1- e-v) Ve-i^ (y{AJ, pV - j) -YiA,j,-l) 



Ca{1 -e-v 



3=0 



T,k^o^k sS pV-no 

(5.31) 
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Now we maximize the second sum in the last hne above over uq, which obviously means 
putting no = 0. The first sum is a geometric series and cancels the prefactor (1 — e~^^^ ). 
As a result, (5.31) is less than 



The key observation now is that, by Proposition 2.4 (a), for given e > there exists m > 
such that \hmo — hoo\ < £ whenever uiq > m. Then 



E 

tiSA'a 



n e-(^W+VVK \h{n,pV) -h^\^eJ2(ll e 



-£{k)nk ^ 



Using the definition of C\, we find that 



(5.32) 



pV 



with 



e"^^' |AiA,p(dx) - UA{dx)\ ^ e + Ca ^ ( V-^ " ^oo)n(A, iV) 

Af=py-m 

^e + SCA E ^o(A,A^), 

Af=pV-m 



(5.33) 



(5.34) 



Note that = ^ yo(A, A^), which suggests that the last term in (5.33) is small as 

the summation from pV — m to pV contains less terms than that giving . To prove 
this, let k be one of the elements of A* closest to 0. By putting j particles into the mode 
k, we find that for any j, N 



Yo{A,N + j) ^ E -ir^ n e"^^^)"" hn^ ^ e^^'^(^) yo(A, A). (5.35) 



k^O 



Thus for any A, 

oo 

^ i?-iyo(A,Ar) e-^'^W = B-'Yo{A,N)- 



j=0 



-ei'k) 



lb(A,A) 
Be{k) 



Inserting into (5.33), we find that 

-Xx 







e \pA^p{dx) - uj>^{dx)\ !^ e + mB^e{k). 



(5.36) 



(5.37) 



As y ^ oo the second term vanishes since e{k) is continuous at 0. Since e is arbitrarily 
small, the left side vanishes indeed in the limit. 

Back to (5.30). For the first term, we can now follow the proof of Theorem A.l in 
[2]. As above, we isolate the contribution of the zero mode and cancel it with the factor 
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1 - e"^/^ . Thus 

e-^^ Mdx) = Ca/ioo (n e-(^(*=)+V^)"'= hn,}j . (5.38) 



Now 



^ e-(^W+VVK/,„^^ =exp(j;iog J] e" . (5.39) 



By Proposition 2.3 (e), the logarithm in the exponential is equal to 

^ e-{e{k)+^)J ^ = ^-e{k)j ^ - - e-^'^^^+y^' e-":' ds. (5.40) 

Using Proposition |2.3| (e) again, the first term in the inner bracket above is equal to 
Sfc^o^^sZ^^o hm and thus it cancels Ca, while the second converges to —pc as 

a Riemann sum for every s € [0, A], when y — > oo. Thus by dominated convergence in s, 
we obtain 

lim / e"^"" ;UA p{dx) = lim / e"^^ z^A(dx) = /loo e"''<= . (5.41) 

Thus by the general theory of Laplace transformations, /ma,p converges to a delta peak of 
strength h^o at pc- The claim of the theorem then follows from (5.26) for p > pc, and this 
completes the proof. □ 

We prove now that the distribution of the random variable n shows typical behaviour. 
To that end we introduce the three sets 

Ae = {n G AAa,7V : -po| < e} 

B^^S = {ne AfA^N : ^ rifc < eV} 

o<\k\<s (542) 

C,,5,M = {^^ G A/a.at : Y nk<eV}. 

fceAMfcl ^ s 

nk>M 



Proposition 5.4. Under the assumptions of Theorem \3.5[ for any density p we have the 
following. 

(a) For any e > 0, lim pA,py(^e) = 1- 

(b) Suppose that pc < oo. For any e > 0, there exists 6e such that PA^pv{B^^sJ > 1 — e 
for V large enough. 

(c) For any e,6 >0, there exists M^^s such that lim PA.pv(Ce,(5,Mj i) = 1- 

The restriction for finite pc in item (b) should not be there — but we cannot prove the 
claim without it. This is the only reason why Theorem 3.5 does not hold when the critical 
density is infinite. 



Proof. The claim (a) follows from Theorem 5.3 



^^_^PA,pv{A) = ^lim^^A,py(l[po-e,po+e](T^)) = / Mpo~e,po+e]is)Spo{s)ds = 1. (5.43) 
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We now get a bound on the probability of a given occupation number. Recall the 



constant B of Corollary 2.5 

PA,N{nk ^ i) 



1 



Y{A,N) 



E 



-e{k)i 



Y{A,N-i) 
Y{A,N) 



(5.44) 



We also have that '^y\a^n)^ ^ (5.22). Then 



EA,N{nk) = ^ PA,N{nk ^ i) ^ 



B' 



i > 1 



(5.45) 



By Markov inequality, 



0<|fc|<<5, 



e Jo<\k\<5, e^W - 1' 



(5.46) 



The integral converges because the critical density (3.8) is finite. It is possible to choose 
V large enough and Se small enough so that PA,n{B'^ < e. 
For the claim (c), we use 



neAfA,N fceA* 



m\ ^ Y{A,N) 

|fc»|><5 



(5.47) 



For given ki, . . . , km, we have 

E n = 



We can bound the last ratio by B. Then 

paacism) ^ E ^ E E 



y(A,iV-^n, 



n 



-s{ki)ni 



ml ^ ^ Yik.m 

m>l A:i,...,fcmeA* ni,...,nm>A/ i=l 



(5.49) 



We bound the ratio of partition functions by -B, see (5.22 ), and we bound one half of eiki) 
by one half of 



en = mill e(A) > 0. 
|fc|>5 



(5.50) 
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Then, since X] rij ^ eA^, we have 
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|fe|>5n>A/ 



\eQep — Be 



4^0 A/ 



-y — 



|fc|><5 



e(k) _ 1 



(5.51) 



We recognise a Riemann sum which is bounded uniformly in V . If M is large enough 
(depending on e and eoj hence on (5), the term in the bracket is positive and everything 
vanishes in the limit V ^ oo. □ 



We are now equipped for the proof of Theorem 3.5, We use Propositions 5.2 and 5.4 
and also Theorem 12.11 



Proof of Theorem \3.5\ From Proposition 5.2 we can split 



0<|A;|<<5 



\k\^5 



(5.52) 



We treat the cases separately. By Proposition 5.4 (a), we can restrict the sum to n G 
with arbitrarily small e. Then for any r]{V) such that r]{V)/V — > 0, 

lim V p{^^pv{n)^Eno{Ni^r,{v)) = ^^"^ yEp^v{Ni {v)), (5-53) 



which is zero by Theorem |2.1[ On the other hand, for the same reasons we have 

lim V pA,py(n)^E„o(iV (v),^y) = ,lim f -Epoy(iV (y) y) 



pV 



s if ^ s ^ po ) 
po if s ^ Po- 



(5.54) 



Next we use Proposition 5.4 (b) to show that the modes < |fc| < (5 contribute a 
vanishing amount. Indeed, we can find arbitrarily small e and b = such that 

^ PKN{n) Y En,{Ni,pv)^e + ^ ^ PA,iv(n) ^ En,{Ni^pv) ^ 2e. 



0<|fc|<<5 



(5.55) 

Therefore we can neglect those modes in (5.52) without changing the result in the Umit 
V ^oo. 

There remain the modes |A:| > 6. By Proposition 5.4 (c) we can restrict the sum over 
n G C^^s,M^ s- because of the definition of C^^s,AU '^^ S^t 



hm ^ 



PA,pv{n) Y En^:{Nr^{V),pv) 

.,5 



^ e. 



(5.56) 



The estimates obtained above prove the second and the third claim of Theorem |3.5 
— and therefore also the first claim, since the fraction of points in microscopic cycles 
is obviously equal to the total density, minus the fraction of points in mesoscopic and 
macroscopic cycles. □ 
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Appendix A. Thermodynamic potentials with periodic boundary conditions 
We clearly have yP*^^(A, N) ^ Y{A, N), so that 

qTip) ^ QAiP)- (A.i) 

It is thus enough to show that liminf converges to q. First we establish some continuity 
property for . 

Lemma A.I. For any A and any rj > 0, we have 



with B the constant of Corollary \K 

Proof. Let N 
conditions is 



p|A| and M = r?|A|. The partition function with periodic boundary 
yP^^(A, N + M)= n e-"^-^^^'^) hn^ 

Yl e-"'==A('^) hr. 



leA/A.jv+A/ fceA* 



(A.2) 



neArA,N A;eA*\{0} 

To get the second line we restricted the sum over occupation numbers to those with 
no ^ M. We know from Corollary 2.5 that hn^-^^M ^ hr^jB and we get the lemma. □ 



Proof of Theorem 3.4 We can derive an expression for the pressure like in the proof 
of Theorem 3.1 namely 



Pa (p) 

The last sum is less than X^fc ^~ 



1 

~W\ 

-eA(fc) 



E 

n > 1 



^fin-On 



n 



fceA* 



■neA(fe) 



(A.3) 



Then converges to the expression for 



p in Theorem 3.1 by dominated convergence. 

We turn to the free energy with periodic boundary conditions. Let us define 

lim inf q^^^'^ip). 



per 



(P) 



(A.4) 



Suppose that there exists p* such that (f^^i^p*^ < q{p*)- It follows from Lemma A.I that 
there exists an interval / close to p* and an r/ > such that 

qrip)<QAip)-V (A.5) 



for all p ^ I. Using (A.I ), we find that the pressure satisfies 



■'?a(#)+'?1/(#) 



(A.6) 



N >0 



Then the infinite volume limit of p^ (p) is larger than the Legendre transform of q — rjXj, 



hence larger than p{p) for some p < 0. This contradicts the first claim of Theorem 3.4 
This shows that q^^ converges pointwise to q. The uniform convergence on compact sets 



follows from (A.I), the uniform convergence of q^, and Lemma A.I 



□ 
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